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Abstract 

We present results from Subaru/FMOS near-infrared (NIR) spectroscopy of 118 star-forming galaxies 
at z ~ 1.5 in the Subaru Deep Field. These galaxies are selected as [Oll]A3727 emitters at z fa 1.47 and 
1.62 from narrow-band imaging. We detect Ha emission line in 115 galaxies, [O III] A5007 emission line in 
45 galaxies, and H/3, [N II] A6584, and [S Il]AA6716,6731 in 13, 16, and 6 galaxies, respectively. Including the 
[O II] emission line, we use the six strong nebular emission lines in the individual and composite rest-frame 
optical spectra to investigate physical conditions of the interstellar medium in star-forming galaxies at 
z ~1.5. We find a tight correlation between Ha and [Oil], which suggests that [Oil] can be a good star 
formation rate (SFR) indicator for galaxies at z ~ 1.5. The line ratios of Ha/[Oll] are consistent with 
those of local galaxies. We also find that [Oil] emitters have strong [O III] emission lines. The [O III] /[Oil] 
ratios are larger than normal star-forming galaxies in the local Universe, suggesting a higher ionization 
parameter. Less massive galaxies have larger [O III] /[Oil] ratios. With evidence that the electron density 
is consistent with local galaxies, the high ionization of galaxies at high redshifts may be attributed to a 
harder radiation field by a young stellar population and/or an increase in the number of ionizing photons 
from each massive star. 
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1. Introduction 

Recent studies indicate that H II regions in star-forming 
galaxies at z fa 1-3 differ from those in local galaxies (e.g., 
Brinchmann et ah, 2008; Liu et ah, 2008; Shirazi et ah, 
2014; Newman et ah, 2014). The measured ionization pa- 
rameter in star-forming regions at high redshifts is found 
to be higher than in normal star-forming galaxies in the 
local Universe (Nakajima et ah, 2013; Richardson et ah, 
2013; Amorm et ah, 2014b; Holden et ah, 2014; Ly et ah, 
2014a; Nakajima & Ouchi, 2014; Shirazi et al., 2014). It 
is also suggested that H II regions of high - 2 galaxies have 
a harder stellar ionization radiation field (Steidel et al., 
2014) and higher electron densities (Hainline et al., 2009; 
Bian et al., 2010; Shirazi et al., 2014; Shimakawa et al., 
2014). 

The SFR of galaxies is known to increase by ~ 30 times 
from the local Universe up to z « 2.5 at a given stellar 
mass (e.g., Daddi et al., 2007; Ly et al., 2011; Whitaker 
et al., 2014), and the number density of active galactic 


nuclei (AGNs) peaks at z fa 2-3 (Hopkins et al., 2007). 
Moreover, studies suggest that the gas fraction in galaxies 
increases up to ~ 50% at 0 ~ 2 from <10% at z ~ 0 (Leroy 
et al., 2008; Tacconiet al., 2010; Geachet al., 2011). Given 
these facts, it is not surprising that AGN and supernova 
feedback, and a larger contribution of young massive stars 
to the radiation field can alter the conditions of the inter- 
stellar matter (ISM) in high-z galaxies. 

Many studies of the metal content in Hll regions of 
galaxies, and its redshift evolution have revealed a stel- 
lar mass-metallicity relation: more massive star-forming 
galaxies are more chemically enriched. This correlation 
evolves toward lower metallicity (at a given stellar mass) 
at higher redshift (e.g., Tremonti et al., 2004; Savaglio 
et al., 2005; Shapley et al., 2005; Erb et al., 2006; Maiolino 
et al., 2008; Liu et al., 2008; Hayashi et al., 2009; Queyrel 
et al., 2009, 2012; Yuan & Kewley, 2009; Richard et al., 
2011; Rigby et al., 2011; Yabe et al., 2012, 2014; Roseboom 
et al., 2012; Stott et al., 2013; Yuan et al., 2013; Henry 
et al., 2013b; Zahid et al., 2013, 2014; Wuyts et al., 2014; 
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Divoy et al., 2014; Ly et al., 2014a, b; de los Reyes et al., 
2015). However, the redshift dependence of the ISM in 
star-forming galaxies can produce one of the major un- 
certainties in metallicity measurements at high redshifts. 
This is because strong nebular emission line diagnostics 
generally used to estimate the oxygen abundance of high- 
2 galaxies are calibrated with metallicity derived by pho- 
toionization models or metallicity based on electron tem- 
perature ( T e ) in the local Universe (e.g., Pagel et al., 1979; 
McGaugh, 1991; Kewley & Dopita, 2002; Pettini & Pagel, 
2004; Nagao et al., 2006; Kewley & Ellison, 2008). Metal 
emission lines, such as [O II] and [O III] in H II regions, arise 
from collisional excitation of ions with electrons. The in- 
tensity ratios of nebular emission lines are sensitive to not 
only metallicity, but also the ionization state and elec- 
tron density of the ISM. Using diagnostics that are cal- 
ibrated with local galaxies requires that we assume that 
the physical state of the gas in high -2 star-forming galax- 
ies is similar to that of local galaxies. However, recent 
studies are raising doubts that this assumption holds at 
high redshifts. 

In this paper, we focus on z ~ 1.5 star-forming galax- 
ies and investigate the physical state of their ISM and 
how it compares to local galaxies using six strong nebular 
emission lines seen in the rest-frame optical wavelength: 
[Oll]A3727, H/3, [O III] AA4959, 5007, Ha, [Nll]A6584, and 
[Sll]AA6716,6731. [Oil] and [Sll] lines are essential to de- 
rive the ionization parameter and electron density of ISM. 
We obtained NIR spectroscopy for 118 star- forming galax- 
ies emitting an [Oll]A3727 line (hereafter [Oil] emitters) 
at 2 ~ 1.5. Our targets are extracted from two samples of 
[Oil] emitters at 2 « 1.47 and 1.62 in the Subaru Deep 
Field (SDF; Kashikawa et al., 2004), which are identified 
by NB921 (A c =91961, AA=132l) or NB973 (A c =97551, 
AA=202l) narrow-band imaging (Ly et al., 2007, 2012a) 
with Suprime-Cam (Miyazaki et al., 2002) on the Subaru 
Telescope. 

While previous studies at z k 1.5-2. 5 have begun to 
increase spectroscopic samples, most of them only have 
observations based on a limited number of nebular emis- 
sion lines. For example, four lines of Ha, H/3, [O III] and 
[Nil] are measured for a few hundred galaxies at z ~1.4- 
1.6 in several deep fields (Stott et al., 2013; Zahid et al., 
2014; Yabe et al., 2014). Steidel et al. (2014) study 179 
star-forming galaxies at 2 ~ 2.3. However, in these stud- 
ies, limited [O II] information is available. In contrast, all 
of our galaxies have the [O II] luminosity measured by the 
narrow-band imaging. There are a few studies with the 
full rest- frame optical spectra covering [Oil], H/3, [O III] , 
Ha, and [Nil]. Most of the samples are not adequate 
to discuss the difference and diversity statistically and/or 
are biased toward massive galaxies or lensed galaxies (at 
most ~30 galaxies at 2 > 1.5 in Masters et al. (2014) and 
Maier et al. (2014a)). Shapley et al. (2014) have used 
a large sample of 118 galaxies, but it is iL-band selected 
(i.e. , mass-selected) sample, and probes 2 ~ 2.3, compared 
with our samples being [O Il]-selected at 2 ~ 1.5 (Ly et al., 
2012a). 

The outline of this paper is as follows. The NIR spectro- 
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scopic observations and the data reduction are described 
in §2. Identification of emission lines in the individual 
spectra and the measurement of emission-line fluxes are 
also described. In §3, we study the physical state of ISM 
in galaxies at 2 ~ 1.5, focusing on measurements of the 
electron density, dust extinction, AGN contribution, ion- 
ization parameter and metal abundance. We then inves- 
tigate the relationship between Ha and [Oil] fluxes for 
individual galaxies, both of which are widely used as star 
formation indicators (e.g., Kennicutt, 1998). In §4, we dis- 
cuss possible mechanisms causing any differences between 
[Oil] emitters at 2 ~ 1.5 and star- forming galaxies in the 
local Universe. We also compare our results with the pre- 
vious studies on star-forming galaxies at similar redshifts. 
In §5, concluding remarks are provided. All magnitudes 
are provided on the AB system (Oke & Gunn, 1983), and 
a cosmology consisting of h = 0.7, il m = 0.3, and Ha = 0.7 
is adopted throughout this paper. 


2. Observations and Data 

2.1. Star-forming galaxies at z ~ 1.5 

2.1.1. Observations 

Our samples in the SDF (13 h 24 m 38f9, +27 d 29 m 25?9) 
contain 933 NB921 [O II] emitters at 2 ~ 1.450-1.485 with 
[Oil] flux > 5.8 xlO -18 erg s -1 cm -2 , as well as, 328 
NB973 [Oil] emitters at z ~ 1.591-1.644 with [Oil] flux 
> 1.7 xlO -17 ergs -1 cm -2 (Ly et al., 2007, 2012a). Here, 
we briefly describe the selection of the [O II] emitters and 
the measurement of [O II] fluxes. We refer readers to Ly 
et al. (2007, 2012a) for details. The emission-line galaxies 
are selected to have a color excess in 2 '— NB (NB921 or 
NB973), where the 2 '— NB colors are corrected with i — z' 
color to properly estimate stellar continuum flux at the 
wavelengths of the NB filters. Then, selections with BR c i' 
and R c i'z' color-color diagrams are adopted to identify the 
[Oil] emitters at 2 ~ 1.47 or 1.62. The [Oil] fluxes (-Fjoii]) 
are measured as follows: 


F [On] 


= ANB 


/NB ~ fz> 

1- (ANB/A 2 O’ 


(1) 


where fx and AA are the flux density and full width at 
half maximum (FWHM) of bandpass “A” , respectively. 
The FWHM of 2 ' is 955A. 

To detect the six major nebular emission lines in rest- 
frame optical for these [Oil] emitters at 2 ~ 1.5, NIR 
spectroscopy was conducted with the Fiber Multi Object 
Spectrograph (FMOS; Kimura et al., 2010) on the Subaru 
Telescope. FMOS is a unique NIR spectrograph with 400 
fibers over a wide field-of-view (FoV), 30 arcmin (diam- 
eter). The observations consist of two programs (S12A- 
026 and S14A-018, PI: M. Hayashi). The first run was 
carried out on 7-9 April 2012 to observe 313 galaxies at 
2~1.5 (211 NB921 and 102 NB973 [Oil] emitters) in high 
resolution mode in “H- long” (1.60-1.80/zm), using three 
fiber configuration setups. We give a higher priority to 
galaxies with larger [Oil] flux (Figure 1). Six G- or Re- 
type stars with H = 15.5-17.3 mag are also included in 
each setup as calibration stars. We use the higher spec- 
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Fig. 1 Top: Histograms of [O II] flux measured from 
narrow-band imaging for all of the NB921 and NB973 
[Oil] emitters at z ~ 1.5 (open) and the emitters targeted 
in our Subaru/FMOS observations (hatched). Bottom: 
The distribution of [O II] flux measured from narrow-band 
imaging as a function of stellar mass for [O II] emitters at 
z ~ 1.5. The NB921 [Oil] emitters at z ss 1.47 are shown 
in blue, while the NB973 [Oil] emitters at z ~ 1.62 are 
shown in red. The circles in pale colors show all of the 
[Oil] emitters at z ~ 1.5. The large circles show the tar- 
geted 313 [Oil] emitters in this study. Among them, the 
filled ones show the 118 [Oil] emitters with at least Ha 
or [O III] detected. The histograms on the right show the 
distribution of [Oil] flux for the targets (open) and the 
spectroscopically confirmed sample (filled). The numbers 
signify the percentage of the confirmed galaxies out of the 
observed ones in each bin of [O II] flux. The histograms on 
the top show the distribution of stellar mass of the targets 
and confirmed galaxies. The galaxies for which the stellar 
masses are poorly estimated are plotted as galaxies with 
stellar mass of 10' M 0 . 


tral resolution mode ( R ~ 2200), because it has higher 
system throughput and is more effective at avoiding con- 
tamination from strong OH sky-lines than the low resolu- 
tion mode ( R ~ 600). The iJ-band spectra enable us to 
observe Ha, [Nll]A6584, and [S ll]AA6716,6731. The cross 
beam switching (CBS) method is adopted, where half of 
the available fibers are allocated to the targets and, at 
the same time, the other half of them observe nearby re- 
gions of blank sky. That is, a pair of fibers (called “A” 
and “B”) are allocated to a galaxy, and ABAB dithering 
is performed for optimal sky subtraction. Observations 
were only conducted during the last 1.5 nights, because 
the weather was cloudy for the first one and half nights 
on 7-9 April. The total integration time was 315, 165, or 
150 minutes for 143, 10, or 160 galaxies, respectively. The 
typical seeing was 0. 8-1.0 arcsec in the R-band, which was 
measured for the bright stars near the center of the FoV. 

The second run was carried out on 11-12 April 2014 to 
observe 125 galaxies at z ~ 1.5 (81 NB921 and 44 NB973 
[O II] emitters) in high-resolution mode in “ J-long” (1.11— 
1.35/nn) using two fiber configuration setups. The J-band 
spectra enable us to observe H/3 and [O III] AA4959, 5007. 
In the second run, we targeted galaxies observed in the 
first run for spectral coverage of all emission lines. We 
gave a higher priority to the galaxies with detections of 
Ha in the iJ-band spectrum. Moreover, since the J and 
H spectra were taken in different observing conditions on 
different dates, we allocate fibers to the same calibration 
stars for accurate flux calibration in all observation runs. 
The CBS method is also adopted in the second run. On 
11 April, almost all of the night was affected by thin cir- 
rus coverage, and the seeing was 0. 7-1.1 arcsec. The sky 
cleared on 12 April, however, the seeing was poor for 
the first half of the night (~1.4-1.7 arcsec). Thus, we 
only obtained optimal data under clear sky condition and 
sub-arcsec seeing (0. 7-1.0 arcsec) for one half-night on 12 
April. The total integration time was dependent on the 
target; 540 (390) minutes for 49 (76) galaxies, respectively. 
Note that the data taken under marginal conditions are 
also used in this study after correcting for the loss of flux. 

Figure 1 shows the observed [Oil] flux of the 313 [Oil] 
emitters targeted in the FMOS observations as a function 
of stellar mass (§2.1.4). This figure shows that although 
our target sample is biased toward bright galaxies with 
-Fjon] ^5x 10 -17 erg s _1 cm -2 , compared with the whole 
sample of [O II] emitters at z ~ 1.5, we have observed star- 
forming galaxies with stellar mass covering about three 
orders of magnitude. 

2.1.2. Data reduction 

The FIBRE-pac data reduction package (Iwamuro 
et ah, 2012) is used to reduce the FMOS data. We run 
the FIBRE-pac in a standard manner. One difference is 
that we first reduce the data per each set of succeeding 
“AB” sequence, and then all of the data are combined af- 
ter the flux calibration while weighting each frame based 
on the flux observed through FMOS fibers. Note that the 
calibration includes the correction for fiber flux loss which 
is estimated from a calibration star. The effects of differ- 
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Fig. 2 The spectra in J- (left) and 77-band (right) taken with FMOS for one of the brightest [O II] emitters. The red 
line shows the spectrum of the galaxy, and blue lines show the Gaussian profile fits (§2.1.3). The blue dotted lines 
show the wavelength of individual emission lines. The black line shows the noise level. The hatched regions show 
the areas that are masked for strong OH sky lines. The 2-D spectra in A and B fibers are also shown on top. Seven 
nebular emission lines, H/3, [O III] AA4959, 5007, Ha, [N II] A6584, and [Sll]AA6716,6731, are seen. 

ential refraction can be important at high airmasses and 
with relatively small fibers (1.2” diameter). However, the 
fiber positioner for FMOS is designed to mitigate these ef- 
fects (Akiyama et al., 2008). All of the observations were 
carried out at airmass less than 2.0 and we re-position the 
fibers every 30 minutes during our observations to account 
for the different zenith distances. Furthermore, we make 
sure that the stellar continuum spectrum is well calibrated 
and there is no discrepancy in the flux between the spec- 
tra separated into the J-long and 77-long coverage for the 
stars used in the flux calibration. 

For extended sources, additional aperture correction is 
required to correct for flux loss. Yabe et al. (2012) have 
shown that the covering fraction by the 1.2” FMOS fiber 
is 0.45 for typical galaxies at z — 1.4. We use the same 
aperture correction in our galaxies. Note that the covering 
fraction is not strongly dependent on the size of the point 
spread function (PSF), if the galaxies are more extended 
than 1.2”. As mentioned above, the flux calibration by 
FIBRE-pac corrects for the flux loss for point-sources. 

The 1.2” aperture includes —90% of the total flux for the 
point sources. Here we assume that the FWHM of the 
PSF is 0.85” and a Gaussian profile for the PSF, since 
the seeing in the FMOS observations are —0.8-0. 9” in 77- 
band. That is, we already corrected for the flux loss by 
-10%. Thus, we made aperture corrections of a factor of 
2 (=0.9/0.45) for both the J- and 77-band spectra. 

Figure 2 shows the J- and 77-band spectra for 
an individual galaxy that has the largest observed 
[O II] emission-line flux and all seven nebular emis- 
sion lines (H/3, [O m]AA4959,5007, Ha, [Nll]A6584, and 
[Sll]AA6716,6731) are clearly detected. 

2.1.3. Line detection 

We fit Gaussian profiles to the spectra to measure the 
redshift and the flux of individual emission lines, and es- 
timate the signal-to-noise (S/N) ratio. The Gaussian fit- 
ting is conducted separately in the J and 77 spectra. In 
each band, major nebular emission lines are simultane- 
ously fitted; H/3 and [O IIl]AA4959,5007 for J-band, and 
Ha, [Nll]AA6548,6584, and [S ll]AA6716,6731 for 77-band. 

Note that the width of the Gaussians fits is fixed to the 



Fig. 3 The distribution of redshifts confirmed by FMOS 
spectroscopy for the 118 [Oil] emitters in the SDF. 
The blue histogram shows the redshift distribution of 79 
NB921 [O II] emitters, while the red histogram shows that 
of 39 NB973 [O II] emitters. The response curves of NB921 
and NB973 narrow-band filters are also shown, where the 
horizontal axis corresponds to the redshift from which 
[O II] emission is shifted into the narrow-band filter. 

same value of [OlIl]A5007 or Ha. In addition, the flux ra- 
tios of the [O III] and [N II] doublets are fixed to one-third. 
Each emission line is well fitted with single Gaussian pro- 
file with a velocity up to —100-200 km s _1 in a, i.e. , 
no broad component in the lino, suggesting that contam- 
ination from Type-1 AGN is negligible. Errors on the 
emission-line fluxes are estimated using the standard de- 
viation of 500 measurements of the flux, where in each 
measurement a set of Gaussians is fitted to the spectrum 
to which the error with normal distribution with ler spec- 
tral noise is added. We regard an emission line with S/N 
greater than 3 as a detection, and visually inspect it for 
confirmation. 

We succeed in spectroscopically confirming 118 [Oil] 
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emitters at z ~ 1.5 by detecting some emission lines 
at more than 3 a, among which we have detected Ha 
([O IIl]A5007) for 115 (45) [Oil] emitters. We also have 
detected H/3, [Nll]A6584, and [Sll]AA6716,6731 in 13, 16, 
and 6 galaxies, respectively. The success rate of line de- 
tection is more than 60% for galaxies with observed [O II] 
flux larger than ~5x 10 -17 erg s -1 cnr 2 (Figure 1). The 
[O III] and Ha emissions are the strongest lines in J- and 
H-band spectra, respectively. Even if only a single emis- 
sion line is detected in the spectrum, it can be identified 
as being [O III] or Ha, because of the additional detection 
of [Oll]A3727 from narrow-band imaging. Thus, we can 
investigate the luminosities of all major nebular emission 
lines in the rest- frame optical for confirmed galaxies; i.e. , 
[Oil], H/3, [O III], Ha, [Nil], and [S II] . The emission- line 
fluxes measured for the 118 confirmed galaxies are pro- 
vided in Table 1. 

Spectroscopic redshifts are computed by taking an av- 
erage of the detected lines. Figure 3 shows the redshift 
distribution for the 118 galaxies. We note that there is 
a discrepancy between the peak of the redshift distribu- 
tion of the [O II] emitters and the response function of the 
narrow-band filters. The iJ-long observations do not cover 
wavelengths )$1.607/im in our observations, and also the 
sensitivity declines sharply at the wavelengths >1.72/mr. 
Thus, the Ha emission line is hard to detect at z (Si -45 and 
z >1.62. This is a likely cause for the discrepancy between 
the redshift distribution and the response of the narrow- 
band filters. Given an emission-line flux measured with 
narrow-band imaging, galaxies located further away from 
the redshift corresponding to the peak of the response 
curve will have larger intrinsic line fluxes. That is, the 
redshift distribution shows that Ha lines in galaxies with 
larger intrinsic [O II] flux are likely to be detected. This 
also suggests that it is important to correct the emission- 
line flux for the filter response function using spectroscopic 
redshift to measure the line flux accurately, even if the ob- 
served flux is measured with narrow-band imaging. 

A fraction of observed galaxies has no emission line de- 
tected, which may be due to several reasons. First, as 
described above, we are likely to miss Ha emission lines 
at z (5 1-45 or z > 1.62. Second, the Ha line strength is 
fainter than the detection limit. This would be the case 
for targets with observed [O II] flux lower than ~5x 10 _1 ' 
erg s^ 1 cm -2 . Indeed, since we regard an emission line at 
more than 3<r as a detection, we find some observed galax- 
ies with marginal detection that are not in excess of 3<r. 
Third, the flux of the line is reduced by the masking of OH 
lines. We note that this is the case for the non-detection 
of Ha or [OlIl]A5007 for some galaxies. Finally, it is 
possible that targeted galaxies are interlopers. However, 
there are only two instances in NB921 [O II] emitters where 
the [O II] AA3726, 3729 doublet is resolved and detected at 
z = 3.37 and 3.72. Thus, we conclude that there is low 
contamination rate in our samples of [O II] emission-line 
galaxies at z ~ 1.5. 


2.1. J- Stellar mass estimates for [Oil] emitters 

Multi-wavelength data from UV to mid-infrared are 
available in the SDF. These data allow us to estimate the 
stellar mass for individual galaxies by fitting the spectral 
energy distribution (SED) by population synthesis mod- 
els. Ly et al. (2012a) performed the fitting to the near- 
UV-to-near-IR SED for the [O II] emitters in this study. 
Thus, we briefly summarize the procedure of their SED 
fitting. 

We used Bruzual & Chariot (2003) population syn- 
thesis models, and assumed the Chabrier (2003) initial 
mass function (IMF), exponentially declining star forma- 
tion histories with r=0.01, 0.1, 1.0, or 10 Gyr, and solar 
metallicity. Note that the Chabrier (2003) IMF is adopted 
throughout this paper, although Salpeter (1955) IMF is 
adopted in Ly et al. (2012a). The assumption in metallic- 
ity is justified given gas-phase metallicity measurements 
of near solar values, as we will demonstrate later (§3.5.2). 
We also note that the differences in stellar synthesis mod- 
els are minimal for different metal abundances. The red- 
dening curve of Calzetti et al. (2000) is assumed and dust 
extinction of Ay=Q. 0-3.0 mag is considered. 

As shown in Figure 1, the [Oil] emitters with emission 
line detected by FMOS have stellar masses of 3 x 10 8 M 0 
to 2 x 10 11 Mq. With a range of three orders of magni- 
tude in stellar mass, we are able to study both dwarf and 
massive star-forming galaxies at z ~ 1.5. 

2.1.5. Stacking FMOS spectra 

Since H/3, [N II] and [S II] are generally weak for typical 
z ~ 1.5 galaxies, only one target has all five major nebular 
lines detected (Figure 2). In many cases, Ha line is the 
only emission line detected, as shown in Tables 1 and 2. 
Therefore, stacking analyses are essential for investigating 
the average physical state of ISM in star-forming galaxies 
at z ~ 1.5. 

Table 2 summarizes subsamples for stacking. Sample- 
1 contains all [Oil] emitters for which J- and iJ-band 
spectra are available. Sample-2 consists of three subsam- 
ples divided by stellar mass. The three subsamples in 
sample-2 are also selected from all [O II] emitters for which 
J- and .ff-band spectra are available. Sample-3 consists 
of four subsamples divided by stellar mass, including all 
[O II] emitters with iJ-band spectra, irrespective of J-band 
data. Sample-3 is used for the analyses involving only Ha, 
[N II] , and [S II] . The subsamples are defined so that the 
fraction of the detected emission line is as similar as possi- 
ble, which avoids a large contribution of a biased number 
of strong emission lines to the stacked spectra. 

To make the stacked spectra, we transform the observed 
spectra into rest-frame based on the spectroscopic red- 
shift. Then we compute the weighted mean of the spectra 
based on the noise at each wavelength with 3er clipping. 
All stacked spectra, shown in Figure 4, detect five nebular 
emission lines, H/3, [O III] , Ha, [Nil], and [S II] . The line 
ratios in the individual stacked spectra are shown in Table 
3. 
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Fig. 4 The stacked spectra for the subsamples summarized in Table 2. The sample-1, -2, and -3 are shown in bottom, 
middle, and top panels, respectively. In each panel, the spectra are shifted vertically for illustration purposes and 
the dotted lines show the wavelength of individual emission lines. The average stellar mass, {log (M*/M q)), of the 
subsamples increases upward by 1.5 (1.1) dex in the sample-2 (3), respectively. As galaxies are more massive, the 
[OlIl]A5007 ([Nll]A6584) line is weaker (stronger) relative to Ha. 


2.2. The SDSS data 

In order to understand the physical state of ISM in 
star-forming galaxies at z ~ 1.5, comparisons with lo- 
cal galaxies are useful. The local spectroscopic cata- 
logs are extracted from the MPA-JHU release for the 
Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) 1 
(Kauffmann et al., 2003a; Brinchmann et ah, 2004; Salim 
et al., 2007; Abazajian et al., 2009). We use 86,238 objects 
that have S/N > 3 for [Oil], H/3, [O III] , Ha, [Nil], and 
[S II] with stellar mass larger than 10 8 MQat redshifts of 
0.04-0.1. These objects are distinguished as star-forming 
galaxies, AGNs, or galaxies with non-negligible contri- 
bution from AGNs (hereafter composite objects) on the 
Baldwin-Phillips-Terlevich (“BPT”) diagram ( [O III] /H/3 
vs. [Nil] /Ha; Baldwin et al. (1981)). The boundary de- 
fined by Kewley et al. (2001) is used to discriminate star- 


forming galaxies and composite objects from AGNs, and 
then the boundary defined by Kauffmann et al. (2003b) 
is used to distinguish between star-forming galaxies and 
composite objects (see also §3.3). The number of star- 
forming galaxies, composite objects, and AGNs are 67764, 
10764, and 7710, respectively. 

3. Physical state of ISM in galaxies at 2 ~ 1.5 

Six major emission lines in rest-frame optical are avail- 
able for star- forming galaxies at z ~ 1.5 in this study, 
allowing us to investigate their physical state of ISM 
in detail. We rely on stacked spectra to study the av- 
erage properties, since most individual spectra are not 
deep enough to detect all the nebular emission lines at 
S/N>3. However, we also investigate individual measure- 
ments where available. 

We first study the electron density of Hll regions in 


i 
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z ~ 1.5 star- forming galaxies. Electron density has an im- 
pact on line strengths. The Balmer lines (e.g., Ha and 
H/3) are emitted from recombination of ionized hydrogen 
gas, while [O III] and [O II] are emitted by collisional exci- 
tation to forbidden transition states followed by radiative 
de-excitation. We then correct the emission lines for dust 
extinction. Correction for dust extinction in the nebular 
emission is important for reliable measurement of intrin- 
sic ratios of nebular emission lines (e.g., Ly et ah, 2012b). 
Next, we use several line ratios to examine possible con- 
tributions of AGN to the emission lines (§3.3), ionization 
parameter (§3.4), and metal abundance (§3.5). Ionization 
level is closely related to the intensity of the radiation 
field in the nebular gas, while metal abundance reflects 
various galaxy evolution processes. Finally, we investi- 
gate the relation between Ha and [Oil] (§3.6), both of 
which are widely used as indicators of star formation in 
galaxies, and the relation between stellar mass and SFR 
estimated from Ha (§3.7). The Ha luminosity directly re- 
flects the number of ionizing photons from massive young 
stars, while the [O II] strength depends in several ways 
on the physical state of ISM, such as oxygen abundance, 
ionization parameter, electron temperature, and electron 
density. However, the use of [O II] has grown in investi- 
gating star formation activity in high - z galaxies, since Ha 
redshifts into NIR (MIR) wavebands for galaxies at z> 0.5 
(z >2.5). 

3.1. Electron density 

The intensity ratios of the [O Il]AA3726,3729 and 
[S Il]AA6716,6731 doublets are sensitive to the electron 
density of ISM. Note that these intensity ratios are in- 
dependent of dust extinction. Although we have the [O II] 
luminosity for our full sample, narrow-band imaging does 
not resolve the [Oil] doublet. Instead, we investigate the 
intensity ratio of the [S II] doublet using the stacked spec- 
tra, and individual spectra for bright galaxies. 

We find that the galaxies at z ~ 1.5 have electron den- 
sities (n e ) of ~ 10 2 cm -3 . Figure 5 shows that there is 
no clear dependence of [Sll]A6716/[Sll]A6731 on stellar 
mass of the galaxy, and suggests that electron density 
in z ~ 1.5 star- forming galaxies is consistent with typ- 
ical estimates for local galaxies, i.e. , n e ~10-10 2 cm -3 
(Osterbrock, 1989; Shirazi et ah, 2014). The electron den- 
sity shown in the right axis of Figure 5 is estimated from 
the intensity ratio of the [S II] doublet under the assump- 
tion of T e = 10 4 K, using the IRAF task temden (Shaw & 
Dufour, 1994). We also find that [S II] A6716/ [S II] A6731 is 
not dependent on [O III] /[Oil], suggesting no dependence 
on the ionization state of ISM (Figure 5). Masters et al. 
(2014) measure an electron density of n e ~100-400 cm -3 
for composite spectra of galaxies at (z) = 1.85. Rigby et al. 
(2011) found an electron density of n e = 252^28 cm -3 from 
the ratio of the [Oil] doublet in a z = 1.70 lensed star- 
forming galaxy. The electron density of our [O II] emitters 
at z ~ 1.5 is consistent with those of galaxies at similar 
redshifts, and with little evolution as a function of red- 
shift. 

Some studies have found higher electron densities of 


~ 10 3 cm' 3 from [S II] doublet measurements for lensed 
galaxies at z ~ 2 (Hainline et al., 2009; Bian et al., 2010). 
Hainline et al. (2009) derive electron densities of 320- 
1600 and 1270-2540 cm -3 for the “Cosmic Horseshoe” 
(z = 2.38) and the “Clone” (z = 2.00), respectively. Bian 
et al. (2010) measured electron densities of 1029/tggg 3 and 
1166^855° cm -3 for two components of J0900+2234, a 
z = 2.03 star-forming galaxy. We find that two individ- 
ual galaxies have [Sll] line ratios showing larger electron 
density of 2.3 xlO 3 or 8.9xl0 2 cm' 3 at 10 4 K. The un- 
certainties in our individual measurements are large, but 
our stacks suggest that such galaxies with denser gas are 
probably not in the majority. 

3.2. Dust extinction 

The difference between the observed and intrinsic ra- 
tios of Ha to H/3, i.e., Balmer decrement, can be used 
to estimate the attenuation of nebular emission by dust. 
Assuming an electron temperature of T e = 10 4 K and elec- 
tron density of n e = 10 2 cm -3 for Case B recombination, 
the intrinsic ratio of Ha/H/3 is 2.86 (Osterbrock, 1989). 
In converting the value of Balmer decrement into extinc- 
tion at 6563A, we assume the reddening curve of Calzetti 
et al. (2000). 

Figure 6 shows the observed ratio of Ha/H/3 as a func- 
tion of stellar mass. Note that no correction is made for 
the absorption of stellar continuum at the Ha and H/3 
wavelengths. The stellar continuum is not detected in in- 
dividual FMOS spectra, and it is hard to estimate the 
influence of the absorption of stellar continuum on the 
Balmer lines reliably with only our data. Thus, we assume 
that the absorption is weak enough for the correction not 
to be required, if any (see also, e.g., Steidel et al., 2014). 

For individual measurements, about half of them have 
observed Ha/H/3 ratios that are consistent with the unred- 
dened Case B value of 2.86 within the uncertainties. 
However, many galaxies have line ratios of Ha/H/3 <2.86. 
The reason for the lower line ratio is unclear, however, it 
is possible that the error is larger than what is shown in 
the figure, due to the uncertainty in positioning of indi- 
vidual fibers allocated to the galaxies. While assuming a 
higher electron temperature T e > 10 4 K can slightly re- 
duce the intrinsic Ha/H/3, the ratio only decreases to 2.74 
for T e = 2 x 10 4 K and n e = 10 2 cm -3 (Osterbrock, 1989). 
The results of individual dust extinction measurements 
with both Ha and H/3 available can be biased towards rel- 
atively low Ha/H/3 ratios of 2.86, since H/3 is easier to 
detect with less dust extinction or intrinsically large H/3 
luminosity. 

On the other hand, the stacked spectra appear to fol- 
low the relation between A(Ha) and stellar mass given 
by Garn & Best (2010), although the extinction in the 
highest mass bin is lower than the relation. As suggested 
by the dual narrow-band survey simultaneously targeting 
Ha and [Oil] for star- forming galaxies at z~ 1.47 (Sobral 
et al., 2012), we also confirm that the relation of Garn & 
Best (2010) is applicable to most star-forming galaxies at 
z ~ 1.5 using the Balmer decrement (see also Dominguez 
et al., 2013; Price et al., 2014). With the stacked spec- 
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Fig. 5 The line ratio of [S II] A6716/[S II] A6731 as a func- 
tion of stellar mass (top) and [0 III] /[Oil] (bottom). The 
right axis shows the electron density n e corresponding to 
the line ratio of [S II] doublet at T e = 10 4 K. The black Filed 
circles show the ratio of the individual galaxies with the 
doublet detected at > 3cr. Red pentagon, blue diamonds, 
and green hexagons show the average ratios derived from 
the stacked spectra of sample-1, sample-2 and sample-3, 
respectively. The contours show the distribution includ- 
ing 68%, 95%, and 99% of star-forming galaxies in the lo- 
cal Universe, where the thicker contour includes the lower 
percentage of the population. 
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Fig. 6 The observed ratio of Ha to H/3 (i.e. , Balmer decre- 
ment), as a function of stellar mass. The right axis shows 
extinction at 6563A corresponding to the line ratio of Ha 
to H/3, where the reddening curve of Calzetti et al. (2000) 
is assumed. The black filled circles show the individual 
measurements with the two lines detected at > 3er, while 
blue squares show the average ratios derived from the 
stacked spectra (sample-2). The dotted line shows the 
intrinsic unreddened Case B value of Ha/H/3=2.86. The 
stacked spectrum with all galaxies indicates an average 
line ratio of (Ha/H/3}=3.57, resulting in (A(Ha))=0.63. 
The open circles show the galaxies with H/3 detected at 
2-3cr, while the arrows show lower limits on the ratio using 
the 2cr upper limit of H/3 for galaxies without H/3 detec- 
tions. The solid curves show the relation between A(Ha) 
and stellar mass in the local Universe and its Ha distri- 
bution given by Garn & Best (2010). 

trum of sample- 1, we find that our [Oil] emitters have an 
average dust extinction of (A(Ha)) = 0.63. 

In this paper, unless otherwise specified, we adopt 
the following procedure to correct for dust extinction of 
nebular emission. First, we assume the dust extinction 
curve given by Calzetti et al. (2000). For galaxies with 
both Ha and H/3 detected at S/N>3, the dust extinc- 
tion is estimated from Balmer decrement. However, if 
(Ha/H/3)obs <2.86, we assume zero dust extinction. For 
galaxies without an observed Ha/H (3 ratio, the Garn & 
Best (2010) relation between A(Ha) and stellar mass is 
used. 

3.3. BPT diagnostics diagram 

The BPT diagram that compares [O III] A5007 /H/3 
against [N Il]A6584/Ha (Baldwin et al., 1981) is frequently 
used to classify galaxies and understand whether massive 
young stars, shocks, or AGN power the ionization of the 
gas. We illustrate in Figure 7 the average BPT line ratios, 
as well as, measurements for individual galaxies if all four 
lines are detected at S/N>3. Here, we apply the bound- 
ary defined by Kewley et al. (2001) to identify AGNs, and 
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Fig. 7 The BPT diagram showing the [OlIl]/H/3 ratios 
versus the [N Il]/Ha ratio. The black circles show the in- 
dividual galaxies with all of the four lines detected, while 
red pentagon and blue diamonds show the average ratios 
derived from the stacked spectra of sample- 1 and sample- 
2, respectively. Open squares show the line ratios for the 
stacked spectra of mass-selected star-forming galaxies at 
z ~ 1.4 (Yabe et ah, 2014). The dashed-dotted and dot- 
ted lines show the relations fitted to z ~ 2.3 star- forming 
galaxies in the Keck Baryonic Structure Survey (KBSS; 
Steidel et ah, 2014) and MOSFIRE Deep Evolution Field 
(MOSDEF; Shapley et ah, 2014) survey, respectively. The 
contours show the distribution including 68%, 95%, and 
99% of star- forming galaxies (blue), composite objects 
(green), and AGNs (red) in the local Universe, respec- 
tively, where the thicker contour includes the lower per- 
centage of the population. The solid lines are the bound- 
aries defined by Kauffmann et al. (2003b) and Kewley 
et al. (2001) to distinguish the populations at low redshift. 
The dashed line is the upper limit of the distribution of 
star-forming galaxies at z= 1.5 that Kewley et al. (2013b) 
proposed. 

use Kauffmann et al. (2003b) ’s boundary to distinguish 
star-forming galaxies from composite objects and AGNs. 
The figure shows that none of our galaxies are located in 
the AGN region. This result is consistent with the low 
velocity dispersion of the emission lines, which indicates 
that our galaxies are not Type-1 AGN (§2.1.3). Since 
the FMOS hirers measure integrated light, there is still 
a possibility that the galaxies host an obscured or weak 
AGN. In such a case, the AGN would not be expected to 
contribute significantly to the emission lines we observe. 

Several recent studies have reported an offset on the 
BPT diagram for star-forming galaxies at z ~ 1.5-2. 5, 
compared with that of local galaxies (e.g., Erb et al., 2006; 
Rigby et al., 2011; Dominguez et al., 2013; Newman et al., 
2014; Shapley et al., 2014; Steidel et al., 2014; Yabe et al., 
2014; Zahid et ah, 2014). Figure 7 shows that our [Oil] 
emitters at * ~ 1.5 also have an offset on the BPT diagram, 


which is consistent with previous studies. From a pho- 
toionization modeling perspective, Kewley et al. (2013b) 
derive a redshift-dependent boundary to separate purely 
star-forming galaxies from galaxies hosting an AGN. The 
revised boundary classifies our [O II] emitters at z ~ 1.5 as 
star-forming galaxies (Figure 7). 

Kewley et al. (2013a) discusses the influence of the ion- 
ization parameter, the electron density, and the shape 
of UV spectrum on the BPT line ratios. Kewley et al. 
(2013a) finds that a higher electron density and a steeper 
UV spectral slope can result in both BPT line ratios being 
larger. A higher ionization parameter leads to an upward 
shift in [O III] /H/3 in the range of [N ll]/Ha that our 2 ~ 1.5 
[O II] emitters span. 

It is an interesting question whether the offset in the 
BPT diagram is in [Nil] /Ha or [O III] /H/3. We attribute 
most of this “BPT offset” to a higher [O III] /H/3 ratio at 
^ ~ 1.5 than in the local Universe. This is because at a 
given stellar mass, our [Oil] emitters have [Nil] /Ha ra- 
tios that are similar to or slightly smaller than that of 
local galaxies (§3.5 and see also Stott et al., 2013; Yabe 
et al., 2014; Zahid et al., 2014). This implies that an up- 
ward shift of the [O III] /H/3 ratio is necessarily required 
to match the distribution of the local galaxies with that 
of the [Oil] emitters at z ~ 1.5, at a given stellar mass. 
Furthermore, the absence of strong cosmic evolution in 
the typical line ratio of the [Sll] doublet suggests that 
electron density is not responsible for the BPT offset, as 
shown in §3.1. One possible explanation is a larger ion- 
ization parameter. In the next section, we investigate the 
ionization parameter for the [O II] emitters at z ~ 1.5, and 
discuss whether this interpretation is correct. 

3-4- Ionization parameter 


It is known that some galaxies at low redshifts, so- 
called “Green Peas” (GPs; Cardamone et al., 2009), have 
extremely high [O III] equivalent width. These green 
peas are rare, low-mass systems with low metallicity 
(Cardamone et al., 2009; Amorfn et al., 2010). Recent 
studies have found high -2 galaxies that possess strong 
[O III] emission (Ly et al., 2007, 2014a, b; van der Wei 
et al., 2011; Nakajima et al., 2013; Nakajima & Ouchi, 
2014; Richardson et al., 2013; Atek et al., 2014; Amorfn 
et al., 2014b; Holden et al., 2014; Shirazi et al., 2014). The 
ratios of nebular emission lines of star-forming galaxies at 
2=1 3 on the BPT diagram imply that typical galaxies at 
high redshifts tend to have stronger [O III] emission, than 
local star-forming galaxies (Figure 7). 

The ionization parameter is defined as 


_ Qho 
^ 'IirR'lnH ’ 


(2) 


where Qh° is the number of ionizing photons above the 
Lyman limit per unit time, R s is the Stromgren radius, 
and nn is the hydrogen density. The ratio of [OlIl]A5007 
to [Oll]A3727 is sensitive to the ionization parameter of 
the ISM (McGaugh, 1991; Kobulnicky & Kewley, 2004; 
Nakajima & Ouchi, 2014). The effective temperature of 
ionizing sources also affects the [O III] /[Oil] ratio, since a 
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Fig. 8 The line ratio of [O III] / [O II] , which is sensitive to ionization parameter, as a function of R 23 (left) and stellar 
mass (right). The black circles show the ratio of the individual galaxies with all of the required lines detected, while 
red pentagon and blue diamonds show the average ratios derived from the stacked spectra of sample- 1 and sample-2, 
respectively. Different populations up to z ~ 3 are also shown; LAEs at z =2-3 (open squares: Fosbury et ah, 2003; 
Erb et al., 2010; Christensen et ah, 2012b; Nakajima et ah, 2013), LBGs at z =2-3 (open triangles: Pettini et ah, 2001; 
Maiolino et ah, 2008; Mannucci et ah, 2009; Richard et ah, 2011; Wuyts et ah, 2012; Belli et ah, 2013), star-forming 
galaxies at 0= 1-2 (open stars: Queyrel et ah, 2009; Rigby et ah, 2011; Christensen et ah, 2012b), star-forming galaxies 
at z <1 (asterisks: Savaglio et ah, 2005; Maier et ah, 2005; Henry et ah, 2013a; Ly et ah, 2014a, b), and GPs (open 
circles: Cardamone et ah, 2009). The contours show the distribution including 68%, 95%, and 99% of star-forming 
galaxies (blue), composite objects (green), and AGNs (red) in the local Universe, respectively, where the thicker 
contour includes the lower percentage of the population. 


harder UV photon can ionize 0 + (Perez-Montero & Diaz, 
2005). 

The left panel of Figure 8 shows the ra- 
tio of [O III] /[Oil] against the R 23 index, 
([Oll]A3727+[Om]AA4959,5007)/H/3, for galaxies with 
the required line detected. In the local Universe, ~ 0.1% 
of SDSS galaxies have [O III]/ [0 11] > 1.0, and a ratio of 
[O III]/ [O II] > 3.0 is seldom seen. Compared with local 
star-forming galaxies, almost all of our [Oil] emitters 
have higher [O III] /[Oil] ratios of > 1 at all R 2 3 values. 
The [O III] /[Oil] ratios of our [Oil] emitters at z ~ 1.5 
are larger than those of typical galaxies at z < 1, and 
consistent with those of galaxies at z = 1-2 (Queyrel 
et al., 2009; Rigby et al., 2011; Christensen et al., 2012b), 
as well as, Lyman a emitters (LAEs) and Lyman break 
galaxies (LBGs) at z ~2-3 (Richardson et al., 2013; 
Amorin et al., 2014b; Nakajima & Ouchi, 2014). 

The right panel of Figure 8 shows the dependence of 
[O III] /[Oil] on stellar mass. The contours and asterisks 
in the figure suggest that there is a similar dependence of 
[O III] / [O II] on stellar mass for galaxies at z = 0.04-0.1 and 
for star- forming galaxies at z < 1 (Henry et al., 2013a; Ly 
et al., 2014a, b). At a given stellar mass, the [O III]/ [O II] 
ratios of our [O II] emitters are larger than those of typ- 
ical galaxies at 2 < 1, and consistent with those of most 
galaxies at 2 = 1-2 as well as LBGs at z ~2-3. Massive 


star-forming galaxies with stellar mass of ~ 10 1 o Mq at 
z ~ 1.5 appear to have [O III] /[Oil] ratios that are similar 
to less massive star-forming galaxies with stellar mass of 
~ 1O 9 M 0 in the local Universe. 

We measure the ionization parameter, q in units 
of cm s^ 1 , from the line intensity ratios of R 23 and 
O32 ([O lll]AA4959, 5007/[Oll]A3727), using the equa- 
tions given in Kobulnicky & Kewley (2004) where the 
photoionization model of Kewley & Dopita (2002) is 
used. [Oil] emitters at z~ 1.5 have (log(g/cm s _1 ))=8.1, 
while star-forming galaxies in the local Universe have 
(log(g/cm s -1 )) ~ 7.3 (e.g., Dopita et al., 2006; Nakajima 
& Ouchi, 2014). Here, we assume that [Oil] emitters fol- 
low the high-metallicity branch of the relation between 
metallicity and R 23 (see §3.5). Even if the low-metallicity 
branch is assumed, average ionization parameter of [O II] 
emitters is (log(g/cm s _1 ))=7.7. Thus, we conclude that 
[O III] / [O II] ratios of [Oil] emitters at z ~ 1.5 may be in- 
creased by higher ionization parameters. We discuss pos- 
sible mechanisms causing the high ionization parameter 
in high-z galaxies in §4. 

3.5. Gas-phase metallicity 

Gas-phase metallicity derived from electron tempera- 
ture measurements are considered to be the most reli- 
able. This is because metals aid in the cooling of gas by 
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optically-thin radiation, and thus electron temperature is 
sensitive to the gas-phase metallicity. However, such mea- 
surements are difficult to obtain, because the emission line 
allowing us to estimate the electron temperature is too 
weak to be detected for individual galaxies 2 . Therefore, 
several “strong-line” emission-line diagnostics, which are 
calibrated with local galaxies, are widely used to derive 
gas-phase metallicity of high -2 galaxies. Here, we will 
assume that these diagnostics are applicable for high -2 
galaxies. 

However, we must be cautious in interpreting high -2 
metallicity studies. In particular, the measured ioniza- 
tion parameter appears to evolve with redshift (§3.4). 
Furthermore, recent studies have examined the nitrogen- 
to-oxygen abundance (N/O), and find enhancement rel- 
ative to local galaxies with the same oxygen abundance 
(Teplitz et ah, 2000; Amorin et ah, 2010; Masters et ah, 
2014; Shapley et ah, 2014). 

3 . 5 . 1 . N/O abundance ratio 

The [Nll]A6584/[0 II] A3727 ratio is correlated with the 
nitrogen-to-oxygen abundance ratio in H II regions (Perez- 
Montero & Diaz, 2005; Perez-Montero & Contini, 2009). 
Since the ionizing potential energies for nitrogen and 
oxygen are similar, the line intensity ratio is insensi- 
tive to ionization conditions (Kewley & Dopita, 2002). 
The relation between the N/O abundance ratio esti- 
mated from the [N II] and [O II] electron temperatures 
and [Nll]A6584/[Oll]A3727 is derived by Perez-Montero 
& Contini (2009). 

Figure 9 shows [Nll]A6584/[0 Il]A3727 as a function of 
R23, that is, the nitrogen-to-oxygen (N/O) ratio as a func- 
tion of the oxygen abundance (O/H). The distribution 
of galaxies in the local Universe is overlaid in this fig- 
ure. At a given oxygen abundance, the [O II] emitters at 
2 ~ 1.5 appear to have higher [Nil] /[Oil] relative to lo- 
cal star-forming galaxies, and the line ratios of the [O II] 
emitters at 2 ~ 1.5 are similar to those of local composite 
objects. While local AGNs (red contours) have stronger 
[N II] emission, compared with star-forming galaxies, the 
contribution of an AGN in our [O II] emitters is negligible. 
We made certain that similar results are obtained even if 
[N II] A6584/[S II] AA6716, 6731 is used as a indicator of the 
N/O abundance ratio (Perez-Montero & Contini, 2009), 
instead of [Nll]A6584/[Oll]A3727. 

Given the result that galaxies at 2 ~ 1.5 have higher ion- 
ization parameter than local galaxies, the R 23 value does 
not necessarily measure metallicity evolution. Higher ion- 
ization parameter leads to lower metallicity at a given R 23 
(Kewley & Dopita, 2002), suggesting that oxygen abun- 
dance of high -2 galaxies is lower than that of local galax- 
ies. Nevertheless, the [O II] emitters at 2 ~ 1.5 are likely to 
have higher N/O abundance ratio than local star- forming 
galaxies. The result that we find is similar to recent find- 

2 Some studies have succeeded in measuring or providing a strong 
constraint on the electron temperature (e.g., Kakazu et ah, 2007; 
Hu et al., 2009; Yuan & Kewley, 2009; Rigby et al., 2011; 
Christensen et al., 2012a; Ly et al., 2014a, b; Amorin et al., 
2014a) 
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Fig. 9 The line ratio of [N II] A6584 / [O II] A3727 as a func- 
tion of R 23 (top) and stellar mass (bottom). The line ra- 
tio is correlated with the nitrogen-to-oxygen ratio (N/O), 
while R23 is an indicator of oxygen abundance, O/H. The 
black circles show the ratio of the individual galaxies with 
all of the required lines detected, while the red pentagon, 
blue diamonds, and green hexagons show the average ra- 
tios derived from the stacked spectra of sample- 1, sample- 
2 and sample-3, respectively. The star shows the average 
ratio for composite spectrum at ( 2 ) = 1.85 (Masters et al., 
2014). The right axis shows the N/O ratio based on cal- 
ibration by Perez-Montero & Contini (2009). The con- 
tours show the distribution including 68%, 95%, and 99% 
of star- forming galaxies (blue), composite objects (green), 
and AGNs (red) in the local Universe, respectively, where 
the thicker contour includes the lower percentage of the 
population. 
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Fig. 10 The line ratio of [Nil] /Ha: as a function of stel- 
lar mass. The vertical axis on the right shows the oxy- 
gen abundance corresponding to the N2 index, using the 
calibration given by Pettini & Pagel (2004). The black 
filled circles show the ratio of the individual galaxies with 
the doublet detected at > 3<r. The red pentagon, blue 
diamonds, and green hexagons show the average ratios 
derived from the stacked spectra of sample- 1, sample-2 
and sample-3, respectively. The solid line and stars show 
the relations given by Yabe et al. (2014) and Zahid et al. 
(2014) for star-forming galaxies at z ~ 1.4 and z ~ 1.6, re- 
spectively. The contours show the distribution including 
68%, 95%, and 99% of star- forming galaxies in the lo- 
cal Universe, where the thicker contour includes the lower 
percentage of the population. 

ings of Amorin et al. (2010) for GPs and Masters et al. 
(2014) for star- forming galaxies at z ~ 1.85. Teplitz et al. 
(2000) also found high nitrogen abundance for MS 1512- 
cB58, which is a gravitationally lensed star burst galaxy 
at z = 2.72. 

Figure 9 also shows that there is no significant differ- 
ence in dependence of N/O abundance ratio on stellar 
mass between the [Oil] emitters at z~ 1.5 and local star- 
forming galaxies. This suggests that at a given stellar 
mass, the [Oil] emitters have similar N/O abundance ra- 
tio to local star-forming galaxies, which is the same as the 
recent finding of Amorfn et al. (2010) for GPs. On the 
other hand, Perez-Montero et al. (2013) found evidence 
for redshift evolution in the relation between N/O abun- 
dance ratio and stellar mass for z < 0.4. In particular, 
they argue that galaxies at z =0.2-0. 4 are offset in N/O 
abundances from local galaxies by -0.14±0.31 dex. While 
this is in disagreement with our result for [O II] emitters 
at 2 ~ 1.5, we note that the selection functions of our sam- 
ple and Perez-Montero et al. (2013) are different, and it 
is not clear if these differences may affect N/O abundance 
measurements. 

Perez-Montero & Contini (2009) found that it is likely 
that some star-forming galaxies are classified as compos- 


ite objects due to high [N Il]/Ha ratio resulting from high 
N/O abundance ratio. The SDSS subsample of star- 
forming galaxies may be biased toward star-forming galax- 
ies with low N/O abundance ratio, which can lead us 
to underestimate the N/O abundance ratio of local star- 
forming galaxies. This suggests that it is possible that 
at a given metallicity the difference in N/O abundance 
ratio between the [O II] emitters at z ~ 1.5 and local star- 
forming galaxies is smaller than what is shown in Figure 9. 
Also, this may result in the slight evolution of the relation 
between N/O abundance ratio and stellar mass, as seen 
by Perez-Montero et al. (2013). However, Figure 9 shows 
that the line ratios of [N II] A65S4 / [O II] A3727 are similar to 
those of local composite objects, although in the BPT dia- 
gram (Figure 7), most [O II] emitters at z ~ 1.5 are located 
in the star-forming region defined by Kauffmann et al. 
(2003b). Thus, we cannot completely rule out the possi- 
bility that the [Oil] emitters at z ~ 1.5 have higher N/O 
abundance ratio at a given metallicity than local star- 
forming galaxies, while having similar N/O abundance 
ratio at a given stellar mass. The trends of N/O abun- 
dance ratio against oxygen abundance and stellar mass, 
as shown in Figure 9, appear to be due to the evolution 
of mass-metallicity relation between the [O II] emitters at 
z ~ 1.5 and local star-forming galaxies. 

3.5.2. Oxygen abundance from [Nll]X658A/Ha 

Figure 10 shows the [N Il]A6584/Ha ratio as a function 
of stellar mass. The N2 index, log([Nn]A6584/Ha), is 
widely used as a proxy for oxygen abundance at high-z 
(e.g., Denicolo et al., 2002; Pettini & Pagel, 2004). Since 
different diagnostics can lead to different metallicity esti- 
mates (Kewley & Ellison, 2008), one of the merits of using 
the N2 index is ease of a fair comparison between our re- 
sult and previous studies (e.g., Erb et al., 2006; Hayashi 
et al., 2009; Stott et al., 2013; Yabe et al., 2014; Zahid 
et al., 2014). Figure 10 shows that there is a dependence 
of the [Nil] /Ha ratio on stellar mass for both individ- 
ual galaxies and the average from stacked spectra. The 
[N Il]/Ha distribution for individual galaxies at z ~ 1.5 is 
similar to that of local galaxies. We note that individ- 
ual measurements are biased because of requiring [N II] 
emission. Thus, the results from stacked spectra (green 
symbols in the figure) show the proper, average relation 
for [O II] emitters at z ~ 1.5. Our [O II] emitters at z ~ 1.5 
show steeper increase of the line ratio with the stellar 
mass, compared to local galaxies. The difference in the 
[N Il]/Ha ratio between the two populations is 0.1 dex at 
the highest stellar mass bin, while the difference is 0.4 dex 
at the smallest stellar mass bin. 

Higher ionization parameter leads to a lower N2 index 
at a given metallicity (Kewley & Dopita, 2002). With ev- 
idence that there is a redshift evolution of the ionization 
parameter, the difference of the N2 index between various 
galaxy populations at different redshifts does not neces- 
sarily imply a metallicity difference. Since less massive 
[Oil] emitters at z ~ 1.5 tend to have higher ionization 
parameter, the metallicity for less massive galaxies may 
be more underestimated. We also find that the nitrogen- 
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Fig. 11 The flux of Ha is shown against [Oil] flux for [Oil] emitters at z ~ 1.5. In the left panel, Ha flux is not 
corrected for dust extinction (i.e. , observed flux), and [Oil] flux is measured from narrow- band imaging (i.e. , not 
corrected for anything). The middle panel is the same as the left one, but for [Oil] flux being corrected for the filter 
response function based on the spectroscopic redshift. In the right panel, Ha and [O II] are both corrected for dust 
extinction (see the text in §3.2 for the details of the procedure of dust correction). 


to-oxygen abundance ratio for [Oil] emitters at z ~ 1.5 
is possibly higher than that of local galaxies, which can 
lead to overestimating the metallicity (Perez-Montero & 
Contini, 2009; Morales-Luis et ah, 2014). Thus, all we can 
conclude from our data is that there is a mass-metallicity 
relation that extends toward M* ~ 10 9 M 0 at z ~ 1.5. 

Photoionization modeling allows us to take into account 
the dependence of line ratios on the ionization parame- 
ter in estimating the metallicity (e.g., McGaugh, 1991; 
Chariot & Longhetti, 2001; Kewley & Dopita, 2002). 
The I ?23 index and the relationship between metallic- 
ity, ionization parameter, and the line intensity ratios, 
which is derived from photoionization models, enable es- 
timates of oxygen abundance without the concerns de- 
scribed above (Kobulnicky & Kewley, 2004), although at 
a given R 23 value there are two metallicity solutions of 
low and high metallicities. We will estimate oxygen abun- 
dance more reliably with R 23 index, and then discuss the 
mass-metallicity relation and the correlation between stel- 
lar mass, metallicity, and SFR for our z ~ 1.5 [Oil] emit- 
ters in a forthcoming paper. 

3.6. The ratio of Ha to [Oil] 

Hayashi et al. (2013) studied the relation between Ha 
and [Oil] for galaxies at z = 1.47, and concluded that 
the luminosity of [O II] is well correlated with that of Ha, 
and thus [O II] can be used to estimate the star formation 
activity even at z = 1.47 (see also Lee et ah, 2012; Sobral 
et al., 2012). However, this study is based on the stacking 
analysis of narrow-band Ha images for [O II] emitters at 
z = 1.47. Therefore, it is worth investigating the relation 
between Ha and [Oil] for 115 individual galaxies at z ~ 
1.5. 

Figure 11 shows the Ha flux as a function of [Oil] flux 
for individual galaxies. The observed fluxes of Ha and 
[Oil] are shown in the left panel. There is a clear cor- 
relation between Ha and [O II] with a low dispersion in 
the Ha/[Oll] ratio of a = 0.17. However, as discussed 
in §2.1.3, the [Oil] observed fluxes require corrections for 
the response function of narrow-band filter. The mid- 


dle panel demonstrates the importance of the correction, 
which ranges from 1.0 to 4.7, depending on the redshift of 
the [Oil] emitters (Figure 3). It shows that the correla- 
tion between the Ha and [O II] observed fluxes is slightly 
tighter, a = 0.15. We also investigate the intrinsic ratio of 
Ha/[Oll] by correcting for dust attenuation (§3.2). The 
right panel shows again the tight correlation where the 
ratio is 1.64 ±0.19, which is consistent with the line ratio 
of local galaxies (Kennicutt, 1998; Moustakas et al., 2006) 
and with z ~ 2 Ha-selected galaxies (Lee et al., 2012). 
Note that the luminosity range of Ha and [O II] covered 
by individual galaxy samples is different by more than 
an order of magnitude between at z ~ 1.5 and at z ~ 0. 
Nevertheless, the tight correlation between Ha and [O II] 
at z ~ 1.5 implies that the [Oil] luminosity can be used 
to estimate star formation activity in galaxies at 2 ~ 1.5, 
similar to that of local galaxies. This is consistent with 
the result of Ly et al. (2012a) who compared [Oil] and 
UV SFRs in z ~ 1.5 [Oil] emitters. 

Even if the ratio of Ha/ [Oil] is corrected for dust ex- 
tinction, there is a non-negligible dispersion (<r = 0.19; 
Figure 11). Since the Ha/ [Oil] ratio can depend on the 
physical conditions in Hll regions (e.g., Moustakas et al., 
2006; Weiner et al., 2007; Gilbank et al., 2010; Maier et al., 
2014b), we compare the ratio to various properties. Figure 
12 shows that the Ha/ [Oil] ratio depends on stellar mass 
and the ionization parameter, [O III] / [O II] (Spearman’s 
rank correlation coefficient is p = — 0.43 and p = 0.61, 
respectively), while there seems to be no strong corre- 
lation between Ha/[Oll] and [Nll]/Ha, i.e., gas metallic- 
ity (p = 0.06). Thus, Figure 12 suggests that less massive 
galaxies or galaxies with higher ionization parameter tend 
to have larger Ha/ [Oil]. 

Given the earlier result that less massive galaxies have 
larger ionization parameter (§3.4), the dependence of 
Ha/ [Oil] on both stellar mass and [O III] /[Oil] appears 
to be related to the stellar mass-ionization parameter de- 
pendence. Since a higher ionization parameter results in 
a larger 0 ++ ionic population relative to 0 + , it is natu- 
ral to expect a weak [O II] emission for galaxies with high 
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Fig. 12 The dust-corrected Ha/[0 II] as a function of stel- 
lar mass (bottom), [Nll]/Ha (middle), and [O III]/ [O II] 
(top). The solid line shows the Ha/[Oll] ratio of 1.64. 
The number in the upper right corner of each panel 
shows Spearman’s rank correlation coefficient. The ra- 
tio of [N Il]/Ha is a proxy of metallicity, while the ratio of 
[O III] /[Oil] is a proxy of ionization parameter. The ra- 
tio of Ha/ [Oil] shows a significant dependence on stellar 
mass and ionization parameter. 

ionization parameter. Indeed, we notice that at a given 
Ha flux, galaxies with larger [O III] /[Oil] ratio tend to 
have lower [O II] flux. This suggests a dependence of the 
Ha/ [Oil] ratio on the ionization parameter. 

The mass-metallicity relation may cause the stellar- 
mass dependence of Ha/ [Oil]. However, the relation be- 
tween metallicity and R 23 index (e.g., Pagel et al., 1979) 
implies that the emission from oxygen ions ( [O II] and 
[O III]) becomes weaker against Balmer emission (H/3) 
with increasing metallicity. We here assume that the 
galaxies follow the high-metallicity branch of the relation 
between metallicity and R 23 index. If this is true, the 
mass-metallicity relation should lead to a positive correla- 
tion between Ha/[0 II] and stellar mass, which is contrary 
to what we observe. 

We conclude that although the intrinsic ratio of 
Ha/ [Oil] for star-forming galaxies at z ~ 1.5 is consis- 
tent with that of local galaxies, the ionization parameter 
has an impact on the Ha/ [Oil] ratio in Hll regions. 

5.7. Correlation between SFR and stellar mass 

It is well known that normal star-forming galaxies at 
each redshift follow a correlation between SFR and stel- 
lar mass from 2 ~ 0 up to z ~ 3, which is sometimes 
called the “main sequence of star formation” or the “star- 
forming sequence” (e.g., Daddi et al., 2007; Noeske et ah, 
2007; Salim et al., 2007; Peng et al., 2010; Whitaker et al., 
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SFR as a function of stellar mass. The SFRs 
are derived from Ha luminosities using the conversion of 
Kennicutt (1998), with corrections for dust attenuation 
as discussed in Section 3.2. The black filled circles show 
the [Oil] emitters at z ~ 1.5. Blue diamonds, and green 
hexagons show the average SFRs derived from the stacked 
spectra of sample-2 and sample-3, respectively. The solid 
line is derived from a linear fit to the individual [O II] 
emitters; log(SFR) = 0.38 x log (M* /1O 1O M 0 ) + 1.45 with 
a dispersion of 0.23 dex. The contours show the distribu- 
tion including 68%, 95%, and 99% of star-forming galaxies 
in the local Universe, where the thicker contour includes 
the lower percentage of the population. The broken line 
shows the relation for star- forming galaxies at z ~ 1.6 ob- 
tained from Subaru/FMOS spectroscopy (Kashino et al., 
2013), while the dotted lines are the relations for star- 
forming galaxies at z =0. 5-1.0, 1.0-1. 5, 1. 5-2.0 and 2.0- 
2.5, respectively (Whitaker et al., 2014). 

2012, 2014). We investigate the SFRs of the [Oil] emit- 
ters at z ~ 1.5 as a function of stellar mass, where the 
SFRs are derived from the dust-corrected Ha luminosi- 
ties using the conversion given in Kennicutt (1998). Note 
that we adjust the conversion factor to that based on 
Chabrier (2003) IMF to estimate SFRs from Ha lumi- 
nosities. Figure 13 shows that there is a correlation be- 
tween the SFRs and stellar masses, suggesting that al- 
most all of the [Oil] emitters are not extreme galaxies 
but normal star-forming galaxies which follow the main 
sequence at their redshift. From fitting a linear function 
to the individual [O II] emitters, the main sequence (MS) 
is expressed as log(SFR) = 0.38 x log(M*/10 10 M e ) + 1.45 
with the dispersion of ctms = 0.23 dex. Such a “main se- 
quence” of star- forming galaxies at z ~ 1.5 is consistent 
with other studies (e.g., Kashino et al., 2013; Whitaker 
et al., 2014). However, the slope of the relation that we 
have obtained for the [O II] emitters is shallower than other 
studies (Figure 13). This seems due to the selection bias 
for the [Oil] emitters with Ha emission detected by our 
NIR spectroscopy. 
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The local star-forming galaxies are also plotted in 
Figure 13. Although the [Oil] emitters have been demon- 
strated to be typical star- forming galaxies at z ~ 1.5 (see 
also Ly et ah, 2012a), the comparison with the local galax- 
ies indicates that at a given stellar mass, star- forming 
galaxies at z ~ 1.5 have SFRs a factor of ~30 times larger 
than typical local star-forming galaxies. We discuss the ef- 
fect of the enhancement of star-forming activity at higher 
redshifts on the ISM of high - 2 star-forming galaxies in 
§4.1. 

4. Discussion 

Unlike typical local galaxies, we have found that star- 
forming galaxies at z ~ 1.5 have [O III] /[Oil] > 1 and a 
large ionization parameter at all stellar masses. We dis- 
cuss the origin of the large ionization parameter for star- 
forming galaxies at z ~ 1.5, and then compare the line 
ratios for our star- forming galaxies at z ~ 1.5 against sev- 
eral NIR spectroscopic studies of star-forming galaxies at 
1. 5-2.5. 

4-1. The cause of large ionization parameter 

Several mechanisms have been considered to explain the 
high ionization state (e.g., Perez-Montero & Diaz, 2005; 
Kewley et ah, 2013a; Nakajima & Ouchi, 2014; Perez- 
Montero, 2014): the shape of the ionizing radiation field 
which depends on an effective temperature of ionizing 
sources and metallicity, the geometry of nebular gas, and 
the hydrogen density. Since we find the electron density 
of Hll regions in galaxies at z ~ 1.5 is similar to that of 
local galaxies (§3.1), it is unlikely that there is a signifi- 
cant difference in hydrogen density between star-forming 
galaxies at z ~ 0 and 1.5. 

First, we consider the geometry of nebular gas. 
Equation (2) states that a smaller Stromgren radius can 
result in a higher ionization parameter. We begin with 
the simple assumption that the nebular gas has a homo- 
geneous distribution around young stars in H II regions. 
In an ionization-bounded scenario, the radius of Hll re- 
gions is the Stromgren sphere radius, which is determined 
by an equilibrium between photoionization and recombi- 
nation. In the density-bounded case, the radius is deter- 
mined by the size of gas clouds, which is smaller than 
the Stromgren radius. Thus, at a given hydrogen den- 
sity, density-bounded Hll regions are more able to ex- 
plain high ionization parameters than ionization-bounded 
regions. We note that the actual distribution of nebular 
gas is likely to be inhomogeneous. If the distribution of 
nebular gas within H II regions is clumpy, then a density- 
bounded scenario is feasible. However, a more detailed 
discussion on the geometry of nebular gas is not possible 
with the existing seeing-limited data. 

Another possible explanation for the high ionization pa- 
rameter is due to the radiation field. Of course, AGNs can 
result in a higher radiation field, and shocks from galactic 
winds can trigger additional excitation of the ionized gas 
(e.g., Rich et ah, 2011; Fogarty et ah, 2012). However, the 
emission-line ratios and line profiles suggest that an AGN 


hard radiation field is negligible (Figures 4 and 7). Since 
lower metallicity, which is found in star-forming galaxies 
at z ~ 1.5, can lead to a harder radiation field, metallic- 
ity may play an important role in determining ionization 
parameter. Indeed, there is a correlation between ion- 
ization parameter and metallicity (Perez-Montero, 2014). 
However, the left panel of Figure 8 indicates that even 
at a given R 23 index (oxygen abundance), star- forming 
galaxies at z ~ 1.5 have larger ionization parameter than 
galaxies at lower redshifts. 

Both SFR and specific SFR (SFR/M*) at a given stellar 
mass increase with redshift up to z ~ 2-3, and less mas- 
sive galaxies have higher specific SFRs. (e.g., Whitaker 
et al., 2012; Speagle et ah, 2014; Ilbert et ah, 2014; Tasca 
et al., 2014, and see also Figure 13). Higher star-formation 
activity corresponds to an increasing number of ionizing 
photons, which could result in a larger ionization param- 
eter at z ~ 1.5 if the geometry and density of gas do not 
change. 

To investigate how the high star-forming activity af- 
fects the ionization state of ISM, we compare the [O II] 
emitters with local star-forming galaxies with comparably 
large SFRs. They are selected as galaxies with log(SFR) > 
0.38 x log(M*/10 10 M Q ) + 0.8 (Figure 14). Only a very 
small fraction of local galaxies meet this criterion. Figure 
14 shows that many local star-forming galaxies with high 
SFRs are able to have [O III] / [O II] comparable to the [O II] 
emitters at z ~ 1.5, suggesting that high star-formation ac- 
tivity is important for the high ionization state. However, 
local star-forming galaxies with high SFRs have a wide 
range of [O III] /[Oil] ratios, so high star formation activ- 
ity does not necessarily result in a high ionization state. 
The middle panel of Figure 14 implies that star-forming 
galaxies with high [O III] /[Oil] ratios > 1 have different 
metallicity at z ~ 0 and 1.5. Moreover, at a given stel- 
lar mass, there are only a few local galaxies that have 
high [O III] /[Oil] ratios comparable to the [Oil] emitters 
at 2 ~ 1.5. 

As we probe higher redshifts, closer to the formation 
epoch of the galaxies, the ages of galaxies should be 
younger. O-type stars (age<l-10Myr) as well as Wolf- 
Rayet stars (age~3-5Myr) may be present. Also, Stanway 
et al. (2014) shows that high [OlIl]/H/3 ratios with high 
ionization parameter can be reproduced for up to ^100 
Myr by taking into account binary stars, which have a 
significant impact on the evolution of massive and Wolf- 
Rayet stars. Therefore, one possibility is that massive 
young stellar populations harden the radiation held of 
high- 2 : galaxies. 

4-2. Comparison with other high-z galaxies 

Stott et al. (2013) studied 193 Ha emitters at z = 0.84 
and 1.47, which were selected from the High-redshift 
Emission Line Survey (HiZELS). The observation was 
conducted in high resolution mode with Subaru/FMOS. 
Among the Ha emitters, 152 galaxies are at z = 1.47. 
Zahid et al. (2014) studied 162 star-forming galaxies at 
z ~ 1.6 in the COSMOS held, which were also observed 
in high resolution mode with Subaru/FMOS. Among the 
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Fig. 14 The left panel is the same as Figure 13 and the middle and right panels are the same as Figures 8, but for the 
SDSS star-forming galaxies with high SFRs comparable to the [Oil] emitters at z ~ 1.5 also plotted. To compare the 
line ratios between star-forming galaxies with similar SFRs at z ~ 0 and z ~ 1.5, SDSS star-forming galaxies meeting 
the criterion (dashed line in the left panel), log(SFR) > 0.38 x log(M*/lO lo M 0 ) + 0.8, are selected (open circles). 


galaxies with Ha detected in iJ-band spectra, 87 galax- 
ies were observed in J-band. Yabe et al. (2014) studied 
343 star- forming galaxies at z ~ 1.4 in the SXDS/UDS 
field, which were observed in low-resolution mode with 
Subaru/FMOS. These studies have detected Ha, H/3, 
[O III] A5007, and [Nll]A6584 lines for individual galaxies 
and stacked spectra, and then investigated the line ra- 
tios and gas-phase metallicity. The number of galaxies in 
these studies is a factor of 2-3 larger than our sample. 
However, our study is the only one with individual inten- 
sity measurements of [Oil] for 118 galaxies at z~ 1.47 and 
1.62. 

In terms of BPT diagnostics, we note that our result 
is consistent with previous high- 2 : studies. Thus, there 
is a consensus that star-forming galaxies at z ss 1.4-1. 6 
are offset from local star-forming galaxies with different 
[O III] /H/3 and [Nil] /Ha line ratios (Figure 7). 

There is a discrepancy between the stellar mass- 
[Nll]/Ha relation that we measure and those of Zahid 
et al. (2014) and Yabe et al. (2014). At a given stel- 
lar mass, star-forming galaxies in our sample show larger 
[N Il]/Ha ratios (Figure 10). Note that we convert all stel- 
lar mass estimates to a common IMF (Chabrier, 2003), 
where applicable. Therefore, the reason for the disagree- 
ment is unclear. One possibility is a selection bias. Our 
galaxies are selected by their [O II] emission, while these 
studies selected by their rest-frame optical continuum 
(A/, -selected). Stott et al. (2013) find a higher [Nil] /Ha 
ratio as a function of stellar mass. They had selected 
galaxies by their Ha emission. While our result is consis- 
tent with that of Stott et al. (2013) at stellar mass > 10 10 
Mq, the [N Il]/Ha ratios of our galaxies are lower at stel- 
lar mass < 10 10 Mq. Since the sample used in Stott et al. 
(2013) includes galaxies at z = 0.84, it may result in higher 
line ratios, due to the evolving mass-metallicity relation. 
It is also worth considering whether the ionization pa- 
rameter and/or N/O abundance are affected by selection 
effects. However, this investigation requires additional ob- 
servations, as previous studies do not have information on 
the [Oil] emission. 


Large NIR spectroscopic surveys are being conducted 
at higher redshifts ( z « 2. 0-2. 6) with Keck/MOSFIRE. 
Steidel et al. (2014) presented initial results of the Keck 
Baryonic Structure Survey which studies 179 star-forming 
galaxies at z ~ 2.3. Shapley et al. (2014) presented early 
observations of the MOSFIRE Deep Evolution Field sur- 
vey, which studies 118 star-forming galaxies at 2 ~ 2.3. 
They also show that there is an offset on the BPT dia- 
gram from local galaxies (Figure 7 and see also e.g., Erb 
et al., 2006), and that the N/O abundance ratio of galax- 
ies at z ~ 2.3 is likely to be different from local galaxies. 
Therefore, strong [O III] emission and high N/O ratio ap- 
pear to be common in high -2 galaxies at 2 ~ 1.5-2. 5. 

5. Conclusions 

We conducted near-infrared spectroscopy with FMOS 
on the Subaru Telescope for 118 [Oil] emission-line galax- 
ies at 2 ~ 1.5 in the Subaru Deep Field, to investigate the 
physical state of the interstellar medium of typical star- 
forming galaxies at high redshift, using rest-frame opti- 
cal nebular emission lines. The galaxy sample consists 
of [Oil] emitters at z «1.47 and 1.62 selected by NB921 
and NB973 narrow-band imaging with Suprime-Cam on 
the Subaru Telescope (Ly et al., 2007, 2012a). Moderate- 
resolution (i? ~ 2200) spectra in two wavelength regions, 
1.11— 1.35/im and 1.60-1.80/im, were obtained. 

We have detected the Ha emission line in i7-band spec- 
tra for 115 galaxies, and [OlIl]A5007 emission line in J- 
band spectra for 45 galaxies. Among these galaxies, H/3, 
[Nll]A6584, [S ll]AA6716,6731 are also detected for 13, 16, 
and 6 galaxies, respectively. Including the [Oll]A3727 
emission line measured by narrow-band imaging, we use 
six major nebular emission lines in the rest-frame optical 
to investigate physical state of the interstellar medium in 
typical star- forming galaxies at 2 ~ 1.5. 

We have found a tight correlation between Ha and 
[Oil], which suggests that [Oil] can be a good indica- 
tor of star formation activity for galaxies at 2 ~ 1.5. The 
line ratios of Ha/ [Oil] are consistent with those of local 
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galaxies. We have also found that the Ha/ [Oil] line ra- 
tio is dependent on ionization parameter and stellar mass. 
The [O II] emitters suffer dust extinction of (A(Ha))=0.63 
mag, which is estimated from the Ha/H/3 Balmer decre- 
ment. 

We have found that [Oil] emitters have strong [O III] 
emission. The [O III] / [O II] ratios are larger than nor- 
mal star-forming galaxies in the local Universe, suggest- 
ing that ionization parameter is enhanced for typical star- 
forming galaxies at z ~ 1.5. Less massive galaxies have 
larger [O III] /[Oil] ratios. With the result that the elec- 
tron densities in our galaxies are consistent with local 
galaxies, the high ionization parameter of galaxies at high 
redshifts may be attributed to a harder radiation field by 
young stellar populations and/or the increasing number 
of ionizing photon from massive stars. 

At a given oxygen abundance, star-forming galaxies 
at z ~ 1.5 may have relatively higher nitrogen-to-oxygen 
abundance ratio, compared to local galaxies. However, 
no clear difference in N/O abundance ratio is seen be- 
tween the galaxies at z ~ 1.5 and z ~ 0 at a given stellar 
mass. These contrasting results suggest that the evolution 
of mass-metallicity relation is responsible for the observed 
trends in N/O abundance ratio at z ~ 1.5. The differences 
in the physical state of the ISM from local galaxies indi- 
cate that direct measurements of metallicity, not using di- 
agnostics with strong emission lines which are calibrated 
with local galaxies, is required to properly study the mass- 
metallicity relation at high redshifts. 
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Table 1. [Oil] emitters in the SDF that are confirmed by Subaru/FMOS spectroscopy. 


ID 

R.A. 

Dec. 

z 1 

^spec 

[O ll]A3727 t 

H ft 

[O m]A5007 t 

Ha^ 

[Nll]A6584* 

[Sll]AA6716,6731* 

NB921.20415 

201.008072 

+27.211218 

1.465 

2.12+0.05 

— 

— 

2.63+0.31 



NB921.20750 

201.080215 

+27.211807 

1.457 

0.72+0.05 



1.14+0.34 



NB921.21503 

200.954636 

+27.213587 

1.465 

0.75+0.05 

— 

— 

1.46+0.31 



NB921.26406 

201.117203 

+27.227135 

1.474 

1.32+0.07 

— 

— 

2.70+0.80 



NB921.28540 

201.023178 

+27.231798 

1.460 

0.75+0.05 

— 

— 

1.63+0.42 



NB921.41594 

201.143753 

+27.269798 

1.464 

1.00+0.04 



1.61+0.28 



NB921.42750 

201.172729 

+27.270859 

1.488 

1.98+0.14 

— 

— 

1.88+0.23 



NB921.42929 

201.018433 

+27.271896 

1.464 

1.83+0.05 


6.16+1.61 

3.85+0.21 



NB921.49590 

201.027985 

+27.292707 

1.472 

2.58+0.06 

— 

— 

3.87+0.45 

1.24+0.35 


NB921.51600 

201.129227 

+27.295822 

1.469 

2.02+0.05 


1.84+0.34 

3.04+0.28 



NB921.52455 

201.111649 

+27.297915 

1.473 

1.62+0.05 


6.18+0.28 

2.32+0.29 



NB921.53008 

201.323746 

+27.299255 

1.484 

1.27+0.08 

— 

— 

1.70+0.10 



NB921.53257 

200.942825 

+27.301119 

1.477 

3.98+0.03 

3.12+0.01 

13.41+0.19 

8.91+0.23 

0.82+0.21 

0.75+0.14, 0.63+0.13 

NB921.53590 

201.292923 

+27.300709 

1.462 

0.62+0.04 



1.63+0.18 



NB921.53725 

200.904846 

+27.303421 

1.476 

4.76+0.11 



3.40+0.54 



NB921.54675 

200.977539 

+27.304607 

1.462 

1.25+0.13 



3.27+0.42 

2.15+0.42 


NB921.57905 

201.268173 

+27.317118 

1.481 

2.92+0.10 


6.60+0.90 

2.97+0.30 



NB921.59540 

200.910477 

+27.321886 

1.465 

1.47+0.05 



1.59+0.22 



NB921.61282 

201.276123 

+27.327734 

1.483 

1.12+0.04 

1.53+0.35 

1.73+0.37 

2.84+0.13 



NB921.62719 

201.008759 

+27.332539 

1.466 

1.42+0.05 



3.97+0.56 



NB921.63744 

200.943481 

+27.337433 

1.481 

5.36+0.16 

2.17+0.31 

6.13+0.81 

6.22+0.36 



NB921.67247 

201.047211 

+27.349859 

1.469 

0.53+0.04 

— 

— 

0.94+0.20 



NB921.68313 

200.997437 

+27.355768 

1.469 

1.95+0.06 



1.98+0.62 



NB921.69065 

201.120758 

+27.358625 

1.473 

9.15+0.09 


10.40+0.12 

13.20+1.22 

3.88+0.43 

1.71+0.61, 1.13+0.36 

NB921.72036 

201.153427 

+27.368980 

1.455 

1.23+0.05 



2.30+0.17 



NB921.73875 

200.918503 

+27.376816 

1.465 

3.95+0.13 



3.75+0.39 



NB921.74041 

201.194672 

+27.377974 

1.456 

1.69+0.05 


9.84+1.60 

4.53+0.27 



NB921.74192 

201.194168 

+27.377691 

1.456 

0.56+0.02 

1.83+0.31 

1.89+0.31 

1.52+0.21 



NB921.76059 

201.091171 

+27.385765 

1.471 

3.37+0.10 



10.16+0.92 

5.11+0.74 


NB921.80256 

200.976547 

+27.400311 

1.465 

0.48+0.03 

— 

— 

1.08+0.32 



NB921.81619 

200.982117 

+27.405230 

1.466 

0.46+0.03 

— 

— 

0.58+0.15 



NB921.84384 

201.035858 

+27.416105 

1.464 

2.80+0.05 

1.21+0.16 

3.92+0.20 

3.47+0.25 



NB921.84670 

201.296234 

+27.416384 

1.476 

6.12+0.25 



6.31+1.07 



NB921.87249 

201.267166 

+27.426624 

1.456 

2.98+0.08 



7.00+0.53 

2.02+0.28 


NB921.87748 

201.132797 

+27.429152 

1.471 

3.21+0.06 


5.80+1.50 

6.47+0.51 

1.02+0.28 

1.11+0.27, 0.99+0.28 

NB92 1.90789 

201.159622 

+27.440001 

1.487 

1.05+0.08 

1.37+0.07 

4.30+0.20 

3.92+0.16 



NB921.90839 

200.900070 

+27.439596 

1.453 

0.73+0.05 



1.63+0.45 




Physical conditions of the ISM at z ~1.5 


Table 1 (cont’d) 


to 

o 


ID 

R.A. 

Dec. 

^spec 

[Oll]A3727 t 

H pt 

[0 111^5007* 

Ha-f 

[Nll]A6584* 

[S ll]AA6716,6731* 

NB921_91027 

201.097717 

+27.441690 

1.472 

2.33+0.14 

— 

— 

1.54+0.46 



NB921.91549 

201.068466 

+27.442301 

1.470 

2.36+0.09 



4.21+1.39 



NB921.95558 

200.978973 

+27.457714 

1.476 

2.66+0.10 



2.24+0.66 



NB921.96152 

200.954590 

+27.462772 

1.477 

2.07+0.08 



2.40+0.33 



NB921.97477 

201.216415 

+27.464403 

1.485 

0.56+0.06 

0.77+0.10 

3.38+0.14 

1.61+0.16 

0.36+0.07 


NB921.101385 

200.912842 

+27.481104 

1.466 

2.19+0.06 



3.34+0.41 



NB921.102638 

201.060425 

+27.486040 

1.472 

1.03+0.05 

— 

— 

2.20+0.41 



NB921.103150 

201.024399 

+27.487307 

1.481 

1.20+0.07 



1.74+0.19 



NB921.103705 

201.003967 

+27.490063 

1.475 

13.88+0.48 


22.12+4.73 

16.96+2.02 

9.29+1.81 

3.21+1.09, 4.60+1.13 

NB921.104655 

201.087936 

+27.493420 

1.464 

1.22+0.03 


6.37+0.15 

3.40+0.29 

0.61+0.15 


NB921.105983 

201.114471 

+27.500454 

1.485 

1.29+0.10 


3.93+0.31 

2.87+0.39 



NB921_107118 

201.019028 

+27.504000 

1.475 

1.54+0.07 



3.38+0.57 



NB921.109072 

200.994446 

+27.512560 

1.473 

5.21+0.18 



5.93+0.98 



NB921_110460 

201.324646 

+27.516947 

1.476 

3.92+0.12 


4.20+0.85 

4.82+0.86 



NB921.110480 

201.070343 

+27.517275 

1.463 

1.61+0.02 

1.59+0.29 

5.24+0.93 

4.44+0.16 

0.27+0.09 


NB921_110936 

201.257339 

+27.518309 

1.477 

0.51+0.03 

1.12+0.22 

1.88+0.27 

1.73+0.15 



NB921_113063 

201.060104 

+27.529142 

1.475 

6.15+0.18 

— 

— 

2.32+0.51 



NB921_114814 

201.170425 

+27.533533 

1.475 

1.55+0.05 



2.37+0.56 



NB921_116328 

201.170258 

+27.538916 

1.476 

0.99+0.04 


3.50+0.18 




NB921_116643 

201.074005 

+27.541149 

1.481 

1.54+0.11 



3.01+0.13 



NB921_119072 

201.065796 

+27.551952 

1.463 

3.40+0.06 



8.64+0.68 

2.03+0.50 


NB921_119235 

201.058258 

+27.550882 

1.473 

0.54+0.04 


2.64+0.26 




NB921_119758 

201.121780 

+27.553385 

1.468 

1.47+0.02 

2.50+0.21 

12.42+0.54 

4.44+0.31 



NB921_121008 

201.245834 

+27.557402 

1.482 

1.38+0.11 


9.59+0.86 

2.63+0.12 



NB921.122309 

201.245865 

+27.564041 

1.461 

2.01+0.08 



3.37+0.49 



NB921.126941 

201.247879 

+27.579729 

1.462 

1.96+0.06 



4.46+0.55 



NB921.127005 

200.998505 

+27.580166 

1.454 

2.29+0.08 



3.40+0.31 

0.58+0.17 


NB921.127149 

201.096802 

+27.580116 

1.478 

0.86+0.07 

— 

— 

1.23+0.40 



NB921.129002 

200.981079 

+27.589703 

1.454 

0.77+0.04 



2.61+0.45 



NB921.132078 

201.049194 

+27.599440 

1.475 

1.94+0.07 



2.76+0.57 



NB921_133042 

201.080048 

+27.603086 

1.479 

0.80+0.07 



1.34+0.36 



NB921.133516 

201.124435 

+27.609943 

1.473 

2.20+0.02 

2.60+0.35 

2.84+0.71 

7.38+0.50 

2.80+0.23 


NB921.133703 

201.124420 

+27.611273 

1.476 

11.95+0.15 


10.61+1.02 

14.98+0.60 

1.91+0.32 

1.71+0.30, 1.31+0.32 

NB921.137733 

201.037720 

+27.619905 

1.482 

0.43+0.06 

— 

— 

0.90+0.25 



NB921_139980 

201.193390 

+27.629015 

1.475 

0.98+0.08 

— 

— 

1.56+0.48 



NB921.140463 

201.129898 

+27.631865 

1.478 

3.48+0.08 


7.68+0.18 

4.64+0.28 


0.59+0.21, 0.65+0.21 

NB921.148013 

201.166000 

+27.658577 

1.476 

5.18+0.09 


5.41+0.38 

4.75+0.50 
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Table 1 (cont’d) 


ID 

R.A. 

Dec. 

■^spec^ 

[0 ll]A3727 t 

H/3* 

[Om]A5007 J 

Ha* 

[Nll]A6584* 

[Sll]AA6716 ! 6731* 

NB921.154909 

201.143539 

+27.681538 

1.475 

0.49+0.04 

— 

— 

0.79+0.21 




NB921.155168 

201.056686 

+27.683475 

1.477 

2.00+0.06 



2.42+0.52 




NB921.155600 

201.251480 

+27.684908 

1.479 

4.60+0.11 



6.97+0.76 




NB921.159865 

201.242142 

+27.698488 

1.477 

1.57+0.08 



1.51+0.47 




NB921.175932 

201.045105 

+27.753534 

1.466 

0.48+0.03 

— 

— 

0.69+0.17 




NB973.39461 

201.021744 

+27.253485 

1.589 

2.92+0.28 



3.35+0.18 




NB973.53162 

201.155777 

+27.293501 

1.589 

7.87+0.59 


12.54+0.83 

4.96+0.56 




NB973.53375 

201.108231 

+27.291931 

1.595 

0.47+0.13 



1.03+0.34 




NB973.53557 

200.968735 

+27.293095 

1.621 

0.78+0.10 


3.34+0.34 

1.26+0.29 




NB973.57024 

201.082626 

+27.301146 

1.589 

0.86+0.23 


1.32+0.23 

1.49+0.31 




NB973.63423 

201.185745 

+27.323143 

1.590 

2.35+0.33 



2.58+0.53 




NB973.63500 

201.184814 

+27.322449 

1.588 

2.21+0.27 

— 

— 

4.97+0.41 




NB973.65046 

201.151443 

+27.328642 

1.599 

2.07+0.25 



4.33+0.57 




NB973.65525 

201.297684 

+27.330811 

1.616 

0.37+0.09 

— 

— 

0.96+0.13 




NB973.69720 

201.184006 

+27.351362 

1.616 

3.05+0.24 


5.59+0.49 

5.02+0.43 




NB973.71489 

200.983398 

+27.361053 

1.590 

4.01+0.38 



4.31+0.77 




NB973.72086 

200.994644 

+27.363752 

1.590 

2.11+0.29 


3.04+0.21 

2.19+0.33 




NB973.72650 

201.210007 

+27.365982 

1.600 

1.83+0.14 



6.13+0.56 




NB973.72934 

201.216019 

+27.368015 

1.599 

5.08+0.37 



5.06+0.49 




NB973.79083 

201.158401 

+27.401751 

1.588 

1.56+0.27 


4.29+0.22 

2.53+0.27 




NB973.81424 

201.315735 

+27.412804 

1.591 

3.79+0.31 


6.83+1.01 

5.71+0.90 




NB973.85096 

200.967087 

+27.428595 

1.599 

0.74+0.15 

— 

— 

1.23+0.23 




NB973.90084 

201.320602 

+27.450769 

1.591 

0.81+0.21 

— 

— 

1.24+0.29 




NB973.90805 

201.290176 

+27.456099 

1.613 

1.93+0.23 



2.54+0.34 

0.96+0.22 



NB973.91278 

200.929718 

+27.457098 

1.597 

1.84+0.20 



2.68+0.34 




NB973.91531 

201.048340 

+27.458298 

1.632 

0.73+0.14 


3.65+0.14 





NB973.104691 

201.186508 

+27.529871 

1.584 

0.92+0.24 


3.51+0.24 

2.39+0.33 




NB973.104719 

200.941956 

+27.531408 

1.617 

1.88+0.06 

1.63+0.56 

3.50+0.80 

4.02+0.21 




NB973.105584 

201.206970 

+27.536465 

1.606 

2.66+0.20 


6.47+0.82 

5.22+0.43 




NB973.108859 

201.048325 

+27.549301 

1.618 

1.84+0.15 


4.71+0.27 

4.15+0.57 




NB973.109063 

200.925018 

+27.550058 

1.616 

2.57+0.22 



5.17+1.25 




NB973.109503 

200.923508 

+27.552755 

1.599 

2.54+0.19 


7.54+0.25 

6.30+0.69 




NB973.115023 

201.158646 

+27.577322 

1.644 

0.45+0.12 

0.76+0.15 

1.60+0.27 

1.86+0.38 




NB973.116135 

201.195541 

+27.582874 

1.589 

1.86+0.31 



3.38+0.73 




NB973.117371 

201.193390 

+27.589464 

1.585 

1.69+0.38 

— 

— 

1.66+0.51 




NB973.118086 

201.349854 

+27.592388 

1.596 

0.67+0.14 


2.35+0.38 

1.08+0.30 




NB973.123696 

200.982178 

+27.620754 

1.592 

3.12+0.28 



4.61+0.85 





to 
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Table 1 (cont’d) 


to 

to 


ID 

R.A. 

Dec. 

z 1 

^spec 

[On]A3727t H/3* 

[Om]A5007 t 

Ha* 

[Nll]A6584* 

[S ll]AA6716,6731* 

NB973.124743 

200.944931 

+27.624931 

1.580 

4.44+0.80 


1.90+0.45 



NB973.128025 

201.226929 

+27.641090 

1.610 

1.77+0.22 


2.16+0.64 



NB973.128071 

201.285156 

+27.641356 

1.589 

1.43+0.26 

— 

2.15+0.45 



NB973.131961 

201.095764 

+27.658642 

1.585 

2.05+0.38 


2.74+0.46 



NB973.133848 

201.094513 

+27.665018 

1.603 

0.58+0.14 


1.40+0.43 



NB973.136464 

201.200027 

+27.675671 

1.604 

2.19+0.20 

3.66+0.38 

4.13+1.18 



NB973.147549 

201.049255 

+27.726721 

1.598 

1.38+0.17 

5.14+0.78 

2.74+0.34 




Note. — Galaxies without J-band spectrum available are marked with — , and non-detections are indicated with 
^Average of the redshifts measured from detected emission lines. 

*The fluxes are measured from narrow-band imaging, and corrected for the filter response function as well as dust extinction. Fluxes 
are in units of 10~ 16 erg s -1 cm -2 . 

*The fluxes of emission line arc all corrected for dust extinction (see the text in §3.2 for the detail). Fluxes are in units of 10 -16 erg 

— 1 —2 
s cm . 
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Table 2. Subsamples for spectral stacking, which are divided based on stellar mass. 



Coverage 

Ntotal 

Subsample 

N 

log (M*/M 0 ) 

H/3 

[O in] 

Ha 

[Nil] 

[S II] 

sample- 1 

J and 77 

89 


89 

9.26 [8.49, 11.2] 

13 

43 

86 

14 

5 




2a 

9 

8.83 [8.49, 9.00] 

2 

7 

7 

1 

0 

sample- 2 

J and H 

89 

2b 

68 

9.62 [9.00, 10.0] 

9 

30 

67 

8 

3 




2c 

12 

10.3 [10.0, 11.2] 

2 

6 

12 

5 

2 




3a 

47 

9.25 [8.49, 9.45] 

— 

— 

44 

3 

1 

sample-3 

H 

113 

3b 

3c 

30 

22 

9.46 [9.45, 9.70] 
9.85 [9.70, 10.0] 

— 

— 

30 

22 

3 

4 

1 

1 




3d 

14 

10.3 [10.0, 11.2] 

— 

— 

14 

5 

2 


Note. - The median stellar masses are shown in each subsamples, and the values in brackets are the 
minimum and maximum of stellar mass in each subsamples. Sample-1 and 2 contain the [Oil] emitters that 
are spectroscopically confirmed and have FMOS spectra in both J- and 17-band. Sample-3 is the same as the 
others, but the galaxies with FMOS spectra in 17-band available are all included. The number shown is that 
of individual emission lines detected in each subsample. Each stacked spectrum is shown in Figure 4. 


Table 3. Intensity ratios of emission lines for the stacked spectra of individual subsamples which are summarized in 
Table 2. All line ratios are corrected for dust extinction, except for Ha/H/3. 



Subsample 

[O m]/[0 II] a 

[O m]/H/3 b 

i?23 C 

Ha/H/3 d 

[N Il]/Ha e 

[N II] / [O II] f 

[S II] doublet® 

sample- 1 


1.80T0.02 

2.97T0.06 

5.61T0.11 

3.53±0.07 

0.16T0.01 

0.28T0.01 

1.47T0.08 


2a 

4.49T0.19 

4.61T0.26 

7.17T0.40 

3.12T0.18 

0.13T0.02 

0.37T0.06 

1.74T0.43 

sample- 2 

2b 

1.77T0.02 

2.98T0.07 

5.65T0.13 

3.51T0.08 

0.14T0.01 

0.23T0.01 

1.54T0.11 


2c 

1.60T0.05 

1.98T0.09 

3.88T0.17 

3.53±0.15 

0.34T0.01 

0.79T0.03 

1.30T0.11 


3a 

— 

— 

— 

— 

0.07T0.01 

0.14T0.02 

1.19T0.14 

sample-3 

3b 

— 



— 

— 

0.13T0.01 

0.23T0.01 

2.75T0.48 

3c 

— 

— 

— 

— 

0.24T0.01 

0.36T0.01 

1.02T0.11 


3d 

— 

— 

— 

— 

0.30T0.01 

0.42T0.01 

1.45T0.13 


a [0 III] A5007 / [O II] A3727 
b [0 III] A5007 /H/3 

c [O II] A3727+ [O III] AA4959, 5007) /H/3 
d The observed ratio of Ha/H/3 
e [Nll]A6584/Ha 
f [N ll]A6584/[0 ll]A3727 
g [S ll]A6716/[S ll]A6731 


